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Heavy hadron production in e + e~ annihilation is studied in the framework of heavy-quark recom- 
bination mechanism. This mechanism predicts a new type of 3-jet event such as bbq or bbq, which 
is power-suppressed relative to the standard bbg event. Furthermore, heavy quark fragmentation 
function also receives a contribution from this mechanism, which is calculable in perturbative QCD 
as long as the momentum fraction z is not too close to unity. We conjecture that the height of the 
peak of fragmentation function for S-wave heavy hadron scales as a 2 a {m) m/AQCD- 



Most perturbative QCD studies have focused on the 
leading power contribution, guided by the factorization 
theorems [l]. Power corrections usually pose a big chal- 
lenge, either because they may not be factorizable, or 
even if so, many more nonperturbative matrix elements 
need to be involved, thus making phenomenological anal- 
ysis more difficult. 

However, there are situations where power corrections 
play an essential role. Deep inelastic scattering process 
at low Q 2 provides such an example 0] . Another classical 
example is the charm production asymmetries observed 
in various fixed-target experiments [3j, where one must 
resort to a new kind of hadronization mechanism other 
than the fragmentation. 

Inspired by the Nonrelativistic QCD factorization 0, 
a new hadronization mechanism, dubbed heavy- quark re- 
combination (HQR) 5], has been developed to explain 
these charm asymmetries. This mechanism simply ex- 
ploits the fact that after a hard scattering, a heavy quark 
may capture a nearby light parton which is produced in 
the hard scattering and happens to carry soft momen- 
tum in its rest frame, subsequently they can material- 
ize into a heavy hadron, plus additional soft hadrons. 
The inclusive production of heavy hadron in this mech- 
anism can be expressed as a product of hard-scattering 
parton cross section and a nonperturbative parameter, 
which characterizes the probability for the heavy quark 
and light parton to combine into a heavy hadron. A typ- 
ical example for the HQR process in hadron collisions is 
qg — >• B + b + X (We denote 6-flavored meson by B). Ac- 
tually, HQR is analogous to an earlier higher-twist model 
for the light hadron production 0, e.g. qg — > 7rq. In con- 
trast, this model involves a light-cone distribution func- 
tion of 7T convoluted with the hard-scattering amplitude, 
thus more complicated than the multiplicative recombi- 
nation factors in HQR. Roughly speaking, it is due to 
simpler dynamics of a heavy hadron. In a heavy-light 
hadron, most of its energy is carried by the heavy quark, 
and the light partons carry a characteristic momentum of 
order Aqcd- However, the momentum of a light hadron 
is widely distributed among its constituent partons. 

The physical idea of HQR is quite general, and not 
only confined in the hadro- and photoproduction of heavy 
hadron. In this Letter, we will apply the HQR to the 



heavy hadron production in e + e~ annihilation. Heavy 
flavor production in the linear collider program is of great 
interest, because it is a perfect arena to test perturba- 
tive QCD, and perhaps more importantly, because some 
related precision electroweak observables like forward- 
backward asymmetry of b quark may be the sensitive 
probes to the new physics pj. To interpret the experi- 
mental results correctly, a comprehensive understanding 
of hadronization of heavy quark is compulsory. 

To simplify the analysis, we will focus on the bot- 
tom meson production from Z° decay. Extending to 
the heavy baryon production, and considering the e + e~ 
center-of-mass energy off the Z° peak, will be straight- 
forward. 

At order a 2 , bottom can be produced through Z° — > 
bbqq. If each quark independently fragments into 
hadrons, then it represents a regular 4-jet event. Nev- 
ertheless, in a small corner of phase space while q is soft 
in the b rest frame, they can be projected onto an inter- 
mediate bq state with definite color and angular momen- 
tum. After a time of order 0(Aq CD ), this composite state 
hadronizes into a B meson plus soft hadrons. Therefore, 
we end up with a jet containing B from the recombi- 
nation, the recoiling b jet and a light quark jet. The 
corresponding Feynman diagrams are shown in Fig. ^ 
Actually they can be obtained by crossing the diagrams 
for the HQR process of heavy meson photoproduction [j| . 
A striking signal of this anomalous 3-jet event is that the 
third jet is initiated by a light quark, instead of by a 
gluon. The inclusive B production rate from HQR can 
be written 

oT[B} = d ^i Z ° b q( n ) + b + q] p[bq{n) -> B] . (1) 

n 

dT is the perturbatively calculable parton cross section, 
and p n are the recombination factors, where n denotes 
the color and angular momentum of bq. We will con- 
centrate on the S-wave states in this work. These p n 
parameters have been recently defined in terms of non- 
perturbative matrix elements ||. An important prop- 
erty of these parameters is that they scale as Aqcd /mi,. 
Since B is produced inclusively, the quantum numbers of 
bq don't necessarily need match those of B. However, it 
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FIG. 1: Diagrams for the HQR process Z° — > bq(n) + b + 
q. Wavely, thick and thin lines represent Z° boson, bottom 
and light quark, respectively. The shaded blob represents the 
hadronization of bq(n) into B meson plus additional hadrons. 
We have suppressed two additional diagrams in which the 
light quark couples to Z° , since they don't contribute to the 
HQR cross section in the leading 1 /m,b expansion. 



is found that all the existing charm meson asymmetries 
data can be fairly accommodated by a single color-singlet 
parameter p\ — 0.15 0. Since it was fit by using only 
the leading order gg — > cc matrix elements as input, we 
should rescale it by a K factor to take into account the 
NLO effects. We choose p\ — 0.3, which can be readily 
converted into p\ ss p\ m c jm,b = 0.1. 

The momenta of bq, b and q are assigned as p, p and 
q. It is convenient to define two independent dimen- 
sionless variables y = 2p ■ q/M z and y — 2p ■ q/M%. 
They are related to the standard energy fraction variables 
z = 2p°/M z , z = Zf/M z by y = 1 - z and y = 1 - z. 
We also define 7 = ml/M z . In the zeroth order of l/m& 
expansion, we can neglect the mass differences between 
B and b. Followingthe prescription for calculating the 
HQR contributions [5| , we obtain the inclusive B and B 
production rate 



1 dT{B] 
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where To is the lowest order decay rate of Z° to bb. The 
6 mass in Tq is neglected for simplicity. The common 
factor p\ in Eq. and © is a manifestation of heavy 
quark spin symmetry [jj. CP invariance indicates the 
cross sections for bq recombination is identical to those for 
bq recombination. These cross sections are highly asym- 
metric in y and y, which is a generic feature of HQR 
processes. Greater enhancement at small y relative to 
y indicates that light quark q prefers comoving with B 
to with the recoiling b. The ellipsis in Eq. J5J and (J3J) 
represent terms that are less singular at small y. Gen- 
erally, they also have different coefficients for the vector 
coupling and axial-vector coupling, signalling a nonzero 
quark mass breaks the chiral symmetry. Full expressions 
of these cross sections will be given somewhere else. 



In a bbg event, perturbative calculation breaks down 
when the outgoing g becomes soft, which corresponds 
to an almost on-shell b propagator. Similarly, when the 
outgoing q becomes soft, perturbative calculation for the 
HQR process will also become invalid. To ensure that the 
HQR parton process is governed by the short-distance 
scale, besides the restriction on the virtuality of b propa- 
gator, the virtuality of the gluon that splits into qq must 
satisfies (q + Pq) 2 ~ 2q ■ pq ^> Aq CD . The typical mo- 
mentum carried by q inside B is of order Aqcd v, where 
v = p/mb is the B velocity. Thus we impose a constraint 
on the invariant mass of B and q 
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Note this constraint is stronger than y, y ^S> Aq CD /M z 
in the bbg case. 

While these three jets are widely separated in angle, 
the HQR cross sections in Eq. © and are suppressed 
by a factor of a s (Mz)AQcr> , mb / M z ~ 10~ 5 relative to 
bbg. This should be regarded as a rather conservative 
estimate, since the HQR cross sections in Eq. (01 and © 
have an extra 1 jy enhancement relative to the bbg case. 
It is important to investigate the discovery potential for 
this anomalous 3-jet event quantitatively. 

Perturbative QCD implies that the gluon jet has some 
characteristic differences from the light quark jet, e.g., 
with higher multiplicity, broader angular distribution and 
softer fragmentation function [Jjj- But it is impossible 
to distinguish them experimentally on an event-by-event 
basis due to large fluctuation, and only possible on a 
statistical basis. OPAL collaboration has carried out a 
comparison study at the Z resonance for 6, uds jets and 
g jets jll|. They selected about 3,000 symmetric 3-jet 
events, in which the most energetic jet with a mean en- 
ergy 42.8 GeV is tagged to contain a b, the other two low 
energy jets carry equal energies about 24.4 GeV. This cor- 
responds to an event topology where the angle between 
the highest energy 6-jet and each of the low energy jets 
is roughly 150°. These samples were assumed to all be 
the bbg events, but now it is clear there must be a small 
fraction of anomalous 3-jet events misidentified as bbg. 

Since only the b quark was tagged experimentally, we 
need sum over all the B flavors emerging from the re- 
combination. We include all the ground-state B mesons, 
but neglect those orbitally excited mesons and b baryons. 
Assuming SU(3) flavor symmetry, we estimate the ratio 
of the yield for HQR 3-jet events to that for the bbg in 
such a topology: 



J2^dT[Z"^B + b + q] 
dT[Z Q -^b + b + g] 



1.6 x 10 



1.0 x 10" 
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a) , 

b) . 



The input parameters are sin 2 9 W = 0.23, a s (Mz) — 0.12, 
Mz = 91.2 GeV and nib = 4.5 GeV. We have considered 
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two possibilities for the HQR contribution. In a), the 
most energetic jet contains B emerging from the recom- 
bination; in b), the highest energy jet is initiated by the 
recoiling b jet. The constraint in Eq. (@J is satisfied in 
both cases. Because the B and q have smaller invariant 
mass (small y) in the latter case, the 3-jet production 
rate is larger than that in the former case. Including both 
possibilities a) and b), we estimate there are roughly 35 
anomalous 3-jet events out of 3,000 samples, seemingly 
not to be statistically significant. We hope that prospec- 
tive Giga-Z experiment with a much larger number of Z° 
samples will confirm their existences. 

Though the HQR cross section is highly suppressed 
for 3-jet, its magnitude becomes much larger when B 
and q lie in the fragmentation region, i.e., with a small 
invariant mass. Therefore, the dominant contribution to 
the cross section comes from a back-to-back two b jet 
configuration. This motivates us to examine if this HQR 
process also contributes to the b fragmentation function. 

Fragmentation functions can be expressed as universal 
nonperturbative QCD matrix elements [l^. However, 
there is a dramatic difference between light hadron and 
heavy quarkonium. While the fragmentation function 
for q, g turning into 7r, K can only be extracted from the 
data, that for a heavy quark into heavy quarkonium, such 
as -D£_ B , can be calculated in perturbative QCD 13]. 

One may expect the fragmentation of heavy quark into 
heavy-light hadron interpolates between these two cases, 
urzcalculable versus fully calculable in perturbation the- 
ory. Indeed, we will show that it is partially calcula- 
ble in perturbative QCD. Some recent work on heavy 
quark fragmentation function using Heavy Quark Effec- 
tive Theory can be found in Ref . |lj] . 

HQR contribution to the fragmentation function can 
be obtained by integrating the double differential cross 
sections and J3J) over y in the limit 7^0. Naively, 
one may think the HQR process doesn't contribute to 
the fragmentation function, due to the overall factor 7 in 
Eq. and However, this explicit suppression can 
be compensated by the enhancement of cross sections in 
the fragmentation region. To see this, note the lower 
integration boundary of y can be well approximated by 
7j//(l — y). Carrying out the integration, changing vari- 
able y to z, we have 

hqr, v _ 32 p\a 2 s {m b ) z (2 - 2 z + z 2 ) 
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hqr , , 32 pj a 2 s (m b ) z (2 - 2 z + 3 z 2 ) 

- — si ir^r — • (6) 

We have also approximated the upper boundary of y to be 
00, introducing an error of order m 2 /M^. The renormal- 
ization scale here must be taken of order mb, instead of 
Mz- Note the absolute normalization of the above func- 
tions is fixed. Also note these expressions don't depend 
on the couplings of Zbb at all, enforced by that fragmen- 



tation functions must be process- independent. In fact, 
due to the "similarity" between HQR and heavy quarko- 
nium production [jj, these functions can be obtained by 
a shortcut. Simply replacing /s c /m c in the D b ^B c (B*) 
derived in 0] by to absorb the singularity from 

vanishing m c , then one can safely take the limit m c — > 
to recover Eq. (J5J and (JBJ. 

The HQR fragmentation functions for both B and B* 
are not away from zero until z becomes large, and finally 
badly diverges as z — ► 1. As a result, the total HQR 
cross section for inclusive B production is linearly diver- 
gent. The breakdown of perturbative calculation as z 
approaches 1 is not hard to imagine. When z ~ 1, b only 
experiences a slight deceleration, and can readily pick 
up a q from the vacuum to hadronize. This should be 
the dominant contribution to the fragmentation process 
in the endpoint region, which is a genuine nonperturba- 
tive problem, unlikely to be tackled in any perturbative 
framework. This is in a sharp contrast to the B c produc- 
tion, where the charm is too heavy to be excited from 
the vacuum, therefore free from this complication. 

We can derive the range of z in which the perturbative 
calculation is valid. Note that y must be large enough 
so that for this given y, the minimum oiy~jy satisfies 
the constraint Thus we must have 



y 



1 -z > 



A 



QCD 



m b 



(7) 



Note this constraint naturally arises from the fundamen- 
tal parameters in QCD and doesn't depend on the exter- 
nal parameter Mz- 

The z distributions in Eq. JSJ and © are much harder 
than the widely used phenomenological Peterson frag- 
mentation function ^| . This may suggest the soft "vac- 
uum" effects are still important, even in the perturba- 
tively valid z region. However, as was advocated in 
Ref. [Hi ], a heavy quark fragmentation function should 
be decomposed into a perturbative part and a nonper- 
turbative part. While the perturbative component sys- 
tematically takes into account the resummation effects 
associated with b production, the nonperturbative one 
characterizes the universal hadronization of b to B. It 
should be kept in mind, however, the nonperturbative 
part alone is not a physical observable and depends on 
how one defines the perturbative part. We will neglect 
this subtlety. A recent extraction of the nonperturba- 
tive part of fragmentation function shows also a harder 
spectrum than Peterson parameterization [l7| . 

The nonperturbative physics starts to take over when 
z ~ 1 — Aqcd/"^6- We may imagine this is a "matching" 
point where perturbative contributions roughly match 
the nonperturbative contributions. It is also natural to 
assume this is where the peak is located 01 . Therefore, 
the magnitude of D — Aqcd /mb) should be of the 

same order as the maximum of the "true" fragmentation 
function. Based on this assumption, a direct observation 
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from Eq. JSJ and JBJ is that the maximum of fragmen- 
tation function for B is 3 times larger than that for B, 
as expected from spin counting. From Eq. JfjJ, we can 
estimate the peak of the fragmentation function of b to 
B*~ is located at z = 0.93, with a height roughly 
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A, 



QCD 



1.5, 



(8) 



where we take the running coupling a s {mj 3 ) = 0.24, 
Aqcd = 300 MeV. This may be an overestimate for B s , 
since the SU (3) flavor symmetry is broken. If we approx- 
imate the "true" B*~ fragmentation function by Peter- 
son function D(z;eb) with e b = 0.006 |l£ |, and take the 
fragmentation probability fo^B— ~ 0.3 Hj, the "true" 
peak is also around z « 0.93 with a height about 1.7, in 
good agreement with the naive estimate in Eq. JSJ). 

Since p\ cx Aqcd/?™&, we expect the maximum of 
the fragmentation function scales as a 2 s (rnb) mb/AQCD- 
Therefore, the ratio of the maximum of D b ^ B *- to that 
of D c ^ D ,o reads (Note B*~ and D*° are in the same 
heavy meson multiplct) 



a 2 s {mb) mb 
(m c ) m c 



1.3. 



(9) 



where m c — 1.5 GeV and a s (m c ) = 0.37 are chosen. 
The dependence of the ratio on Aqcd is merely from 
the running strong coupling. With e c = 0.05 [l^] and 
fc~>D'° ~ 0-24 19], the Peterson function D c _ D *o peaks 
around z = 0.8, with a height about 0.64. The ratio 
of the maxima of the "true" fragmentation functions is 
1.7/0.64 « 2.6, twice as large as expected from the scal- 
ing law in Eq. 10. Regarding the potentially large l/m c 
corrections, we are content with this crude agreement. 

In summary, we have investigated heavy meson pro- 
duction at e + e~ collider from a new hadronization mech- 
anism. This mechanism predicts a new type of 3-jet 
event, power-suppressed relative to the ordinary bbg jet. 
Even though its discovery is hampered by the experi- 
mental difficulty to distinguish the quark and gluon jets, 
it is possible that future Giga-Z experiment will provide 
overwhelming evidences for this mechanism. 

It is shown that HQR also has a contribution not sup- 
pressed by inverse powers of Mz, but only by Aqcd/to6. 
We identify this with a contribution to heavy quark frag- 
mentation function. We have calculated it in the zeroth 
order of heavy quark expansion. It will be very inter- 
esting to incorporate the l/rrib corrections. The pertur- 
bative calculation breaks down when the energy fraction 
z is close to 1 — AqcD/mi,. Finally we argue that the 
maximum of the fragmentation function for S-wave heavy 
meson may respect a simple scaling law. 

One important issue we haven't addressed so far is how 
to extract the finite power corrections from the power 



divergent total HQR cross section. To do so, we must 
identify and then subtract off all the leading power (frag- 
mentation) contributions in the infrared region. This is 
a nontrivial task and may need invoke a more general 
factorization theorem l20l. 
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